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FOREWORD 


This work, performed at the NASA Langley Research Center, represents 
progress on an ongoing research effort for the development of an analytical model 
for thermoset prepregging processes in autoclave fabrication applications. The 
work was supported under NASA Contract NAS1-16000. Mr. R. M. Baucom (PMB) is 
Technical Monitor for this project. This work was performed by Dr. T. H. Hou of 
Kentron International, Inc. 


ABSTRACT 


A new flow model which describes the process of resin consolidation during 
prepreg lamination has been developed. A parametric study has also been carried 
out to explore the Salient features of model predictions. It is assumed that 
resin flows in all directions are originated from squeezing action between two 
approaching adjacent fiber/fabric layers. In the horizontal direction, a 
squeezing flow between two nonporous parallel plates is analyzed, while in the 
vertical direction a poiseuille type pressure flow through porous media is 
assumed. Proper force and mass balances have been established and solved for the 
whole system which is composed of these two types of flow. 

A flow parameter, Cp, has been defined and shown to be a measure of 
processibility for the curing resin. For a given external load F the responses 
of resin flow during prepreg lamination, as measured by Cp» are categorized 
into three regions: the low Cp region where resin flows are inhibited by the 

high chemoviscosity during initial curing stages; the median Cp region where 
resin flows are properly controllable; and the high Cp region where resin flows 
cease due to fiber/fabric compression effects. Resin losses in both directions 
can be calculated as well. 

Potential uses of this model, including quality control of incoming prepreg 
material, are discussed. Experimental measurements needed to confirm the prepreg 
characterization methods suggested by the model have been planned for future 
research. The parallel -pi ate plastometer will be an appropriate device to carry 
out the experiments. 
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LIST OF SYMBOLS 


AE n 

AE k 

F 

h(t) 

h„ 


Area (cm 2 ) perpendicular to the z-di recti on 
Flow Parameter (cm 2 /dyne) defined in Eq. (11) 

= dp/dz (psi/cm) 

Viscous flow Activation Energy (Kcal/mole) 

Viscosity cure Activation Energy (Kcal/mole) 

External loading (dynes) 

Separation (cm) between parallel plates 
Initial separation (cm) between parallel plates 


1 - h(t)/h A measure of plate separation (Fig. 2) at any instant t 
fraction of the original distance h Q 


h(t) 

K 

k 

k 

cc 

p 

Pa 

P G 

Qi »Q2 

Q r .Q z 

Qt 

R 

r,0 ,z 


= dh/dt (cm/sec) 

Permeability (cm 2 ) of porous material 
Viscosity rate constant 
Material constant (min -1 ) 

Pressure (psi) generated by squeezing action between two 
approaching plates 

Ambient pressure (psi) 

Pressure (psi) absorbed by glass fabric or fiber bundles 
Defined by Eq. (10) 

Volumetric resin flow (cm 3 ) in r and z direction, 
respectively 

Total volumetric resin flow (cm ) 

Universal gas constant (Kcal/mole °K) 

Cylinderical coordinate system 
Curing temperature (°K) 


as a 
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Curing time (sec) 

Velocity (cm/sec) averaged over a small region of space in 
porous material 

Velocity (cm/sec) in r and z direction, respectively 

Characteristic thickness (cm) of porous material where resi 
flows 

Chemoviscosity (poise) 

Initial viscosity (poise) at t = 0 

Material constant (poise) 

Density (gm/cm 3 ) averaged over a region in porous material 



I. INTRODUCTION 


High quality fiber reinforced composites, such as those used in aerospace 
and industrial applications, are commonly processed in autoclaves. During 
processing, the composite materials are subjected to prescribed elevated tempera- 
tures and pressures. Selection of a cure cycle (temperature profile) will 
dictate the kinetics of polymer chemical reactions and the viscosity growth 
profile. In order to maintain control over these parameters an adequate 
chemoviscosity model is required. During processing, increased molding pressure 
results in resin flow both perpendicular and parallel to the fibers. The applied 
pressure helps to consolidate the composite laminates and to squeeze out excess 
resin. An adequate resin flow model for the entire system (laminate/bleeder/ 
breather) is needed to provide a description of the time-dependent laminate 
consolidation process by analyzing the flow pattern, loss of resin, and pressure 
profiles during processing. 

Considerable work has been conducted in the past by numerous researchers in 
searching for a relationship between chemoviscosity and cure kinetics . 1 However, 
the mechanics governing the flow of resin associated with composite laminate 
processing has received little attention. Bartlett 2 developed a theoretical flow 
model for glass-reinforced resin during lamination for multilayer printed circuit 
boards in the electronics industry. His model assumed that the glass fabric 
behaves as a porous media; that the flow over and through the glass bundles 
offers the greatest resistance to flow (i.e., resin flow in perpendicular 
direction is negligible); and that the lamination process is isothermal. A 
poiseuille type of flow was assumed for the horizontal directions. The unique 
features of Bartlett's model are that (i) the specific permeability of the porous 
glass reinforcement is considered a function of fabric geometry; and (ii) the 
effect of layers of glass fabric coming into contact with one another as the 
resin is squeezed out has been taken into account. The model was later compared 
with experimental results by Bloechle 3 *** using a paral lelplate plastometer, 
and was successfully applied to quality control of incoming epoxy B-stage 
prepregs in a manufacturing environment. 

Springer and Loos 5-8 considered resin flows perpendicular and parallel to 
the planes of the composite. Darcy's law for flow through porous materials was 
assumed in the vertical direction, while poiseuille type pressure flow was 



analyzed in the parallel direction. Macroscopic force and mass balances were 
performed in each direction from which the approaching speeds dh/dt for two 
adjacent layers of porous plates separated by a distance of 2h (see fig. 2), and 
the amount of resin loss was derived. Experiments were conducted with Hercules 
AS/3501-6 graphite epoxy prepreg and the model predictions were shown to compare 
favorably with the experimental data for multiple laminates. 

Lindt 9 studied the flow mechanics associated with the formation of fiber 
reinforced composite laminates. The overall flow pattern in the laminate was 
conceptually compartmentalized into units of a flow cell. Each individual cell 
is composed of two types of flow - a squeezing flow in the horizontal direction 
between two fiber layers, and a poiseuille type pressure flow in the vertical 
direction. Both are originated from the squeezing action between two adjacent 
fiber layers. Lubrication approximation is employed to formulate the equations 
of equilibrium. Numerical solutions reveal reasonable patterns of fiber 
distribution, loss of resin, flow field, and pressure profiles in the laminate. 

No experimental comparisons have been reported, however. 

In this report, a newly developed resin flow model for composite processing 
will be described. It is assumed that flows in both vertical and horizontal 
directions originate from the squeezing action between adjacent fiber layers. A 
parametric study has been performed to illustrate the salient features of the 
model predictions with regard to the relationships between process conditions and 
material responses. Only the case of isothermal cure has been considered in this 
study . 
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II. THEORY 


Physical Model 

The schematic diagram of a composite laminate is shown in Fig. 1. A stack 
of 5 layers of prepreg tapes is confined between two steel plates and layers of 
porous materials. The porous materials include those commonly called bleeder and 
breather materials. When a force, F, is applied, the resin flow pattern observed 
by Springer 7 indicates that the first (top) layer moves toward the second one, 
while the resin is squeezed out from the space between the two layers. 
Simultaneously, a similar situation holds for the bottom layer. When the first 
layer reaches the second one, the two layers will move together toward the third 
layer, squeezing resin out of the space between the second and third layers. 

This sequence is repeated for the subsequent layers. One important observation 
noted by Springer is that upon application of the force F, the first (top) plate 
moves downward (or upward for the bottom plate in Fig. 1). A plate further down 
starts moving only when the plates above reach it. 

A typical flow channel between two layers of prepreg tapes is shown 
schematically in Fig. 2(a). Resin is confined between parallel porous plates 
separated by a distance, h. Upon application of a force, F, the resin is 
squeezed outward horizontally and also vertically through the porous media. It 
is assumed that these two flow directions can be decoupled conceptually as shown 
in Fig. 2(b) and (c). Fig. 2(b) illustrates a squeezing flow between two 
non-porous plates separated by a same distance, h, as in Fig. 2(a). A vertical 
pressure flow through porous materials is illustrated in Fig. 2(c) (assuming no 
horizontal flow), and that flow through the porous materials is driven by a 
pressure drop (p - p a ) across a characteristic distance, Z 0 , where p is the 
pressure generated by the squeezing action between plates, and p a is the 
ambient pressure. The same external force, F, is applied in Fig. 2(b) and (c). 

Mathematical Formulation 


A. Vertical Flow Through Porous Media 

For the flow of a fluid through a porous medium, the equation of motion 
can be replaced by Darcy's law 10 



v 0 = . £ (vp - pa). 


(i) 


where the underlined quantities denote vectors. K is the permeability of the 
porous medium, n is the viscosity, V Q is a superficial velocity averaged over a 
small region of space, and p and p are density and pressure, respectively, 
averaged over a region available to flow that is large with respect to the pore 
size. 

For an incompressible liquid and constant K and n, Eq. (1) together with 
the equation of continuity can be reduced to 

V 2 p = 0 (2) 

As a first order approximation, we assume that Eq. (2) is applicable to our 
system (Fig. 2(c)) in a unidirectional flow. Eqs. (1) and (2) can then be 
combined to give 



(3) 


where C z = dp/dz is a constant. The volumetric resin flow in the (vertical) 
z-di recti on is therefore equal to 


Q 7 = AC K / -±- dt 
z z z 0 n(t) 


(4) 


where A z is area perpendicular to the z-direction and n(t) is the 
chemoviscosity of the reacting resin. 

B. Horizontal Squeezing Flow 

We now consider the squeezing flow between two non-porous plates, as 
shown in Fig. 2(b). A cylindrical coordinate system (r,0,z) is chosen for 
convenience. The velocity and pressure fields are assumed V r = V r (r,z), 

Vz = V z (z) and p = p(r) only. The equation of continuity becomes 
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where h = dh/dt is a measure of the speed of two approaching parallel plates 
upon the application of force F, and V z = h at z = h. 

For the velocity field assumed, Bird et al. 11 have shown that the 
r-component equation of motion can be integrated to give 

t 1 - <F> 2 ] < 6 > 

Substituting Eq. (6) into Eq. (5), we have 


dp 3_ n_r 
dr = 2 ,3 


(h + 


V’ 


which can then be integrated to obtain the pressure distribution as 


(7) 


p(r)-p a 



<H * V 0 ) [1 - (f) 2 ] 


( 8 ) 


Eq. (8) represents the squeezing generated pressure drop between two parallel 
plates approaching with a speed h(t). By equating a balance in force, we have 


F 



[(P(r)-p a ) + C z Z Q + P G ] 2irrdr 


(9) 


= *R 2 (C z Z 0 + P G ) -f *4 [nh + KC Z ], 

h 


where F is the applied external force to the plates, and Pq is the pressure 
absorbed by the glass fabric or fiber bundles. Physically it is noted that when 



layers of fibers come into contact with one another as the resin is squeezed 
outward, they begin to carry a portion of the applied load. The average pressure 
applied to the resin is therefore the difference between the average applied 
pressure and the pressure Pq carried by the fibers. Pq was shown to be a 
function of distance h by Bartlett 2 . Eq. (9) can be rearranged to become: 


t 


I 


o 


1 

TtftT 



l 


KC, + ( 


Qz 


)h : 


dh. 


( 10 ) 


with 


Qi = ttR 2 ( c z z o + F q) » 


and 


O2 = ■g *R 4 

For a given fabric subjected to a known load F, Eq. (10) then relates the plates 
approaching characteristics h(t) to the change of chemoviscosity (rate of cure) 
of thermosetting resin. A quantity called the flow parameter, Cp» can 
therefore be defined as: 


C F =/t „W dt ‘ 


(ID 


The significance of Cp will be discussed in the following section. 


Substituting Eq. (7) into Eq. (6), we have 


v r (r.i) - -!£ <h ♦ V [1 - (f> 2 L 


( 12 ) 


6 



and Qp» the volumetric edge resin loss, can be calculated. A more direct 
method is to calculate the total resin loss, Qj, as 

t • 

Q T = - A / h dt, (13) 

o 

and then Q r can be obtained from 

Q r = Q T - Q z » < 14 ) 

where Q z is calculated by Eq. (4). 
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III. Numerical Results and Discussions 


For the purposes of illustrating unique features of model predictions 
discussed previously, a 1080 B-stage glass fabric material used by Bloechle 3 was 
chosen. Permeability, K, of the material is 5.8 x 10“ 12 cm 2 . The 
load-deformation characteristics of the same material has been reported by 
Bartlett 2 and is also reproduced in Fig. 3. 

From Eqs. (10) and (11), we have 


C f = ! o 


dt = - / 


1 

h n KC 7 + (^-jr 2 ^ 3 
0 z Q 2 


dh 


(15) 


For a pair of parallel plates with a radius R = 7.45 cm separated initially by 
a distance h 0 = 9.144 x 10" 3 cm, Eq. (16) can be solved for fixed value of 
pressure gradient, C z , in vertical flow through a porous media as shown in 
Fig. 4, or for a fixed value of external load F, as shown in Fig. 5. The 
quantity l-h/h 0 in the vertical axis represents a measure of processibility or 
resin flow, and is shown to be related directly to the flow parameter Cp. 

Similar characteristics are noted to exist for each curve in Fig. 4. In the 
region of low Cp where values of l-h/h 0 are relatively unchanged, the flow 
characteristic is dominated by chemoviscosity changes in the resin during the 
initial curing stages (see Eq. (15)). In the middle region of each curve where 
values of l-h/h 0 are most sensitive to Cp, resin flows can be properly 
controlled by adjusting the values F; namely, more excess resin can be squeezed 
out for higher applied loads. As the values of Cp increase, a plateau region 
for each curve becomes apparent. In this plateau region, values of l-h/h 0 
become insensitive to changes in Cp again since adjacent fiber or fabric 
layers nest together and are compressed against each other. At this stage fibers 
carry the major portion of the load and the resins cease to flow. 


It can be noted from Fig. 4 that reductions in the external applied load F 
shift the curves to the higher Cp region. The maximum attainable resin flow 
l-h/h 0 is also reduced. Both of these behaviors are consistent with the above 
arguments. Upon application of two different loads with F J <F 2 , the 
corresponding horizontal shift between two curves in Fig. 4 can be calculated by 
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h 3 (F 2 -F 1 )/Q 2 


r 1 . r 2 - f 1 dt = f ' 2 1,,vz 

C F C F / ti nftT ” \ [K C z t (£ljgl,h*] [KC 2 ♦ (Jiljgljh*] 


dh, (16) 


and is shown to be a function of distance, h. 

Fig. 5 shows the effects of the magnitude of C z Z 0 . C 2 Z 0 is the 
pressure drop which drives the resin to flow vertically through the porous 
media. Increasing the magnitude of C z Z 0 has a similar effect to decreasing 
the external load F as discussed in Fig. 4. The horizontal shift between two 
curves of different C z Z 0 values can be calculated by 


t 2 1 

Cp 1 - C F 2 = / — dt 

F F J . n 


h 

/ [ 

h 


0 KCi + (^2l)h 3 

Z 1J2 


KC i + 


-] dh. 


(17) 


for (C z Z 0 ) 1 > (C z Z 0 ) 2 . The amount of shift (Cp 1 - Cp 2 ) due to 

changes in magnitude of C z Z 0 is, however, smaller than that due to changes 

in the magnitude of external load F, as can be seen by comparing Figs. 4 and 5. 

The B-stage epoxy viscosity can be represented by 2 


n(t) = n 0 exp [kt]. 


(18) 


with 


n o = n « exp [ AE n /RT ]» 


and 


k = k a exp [-AE k /RT], 


9 



where n and n 0 are curing and initial viscosity, respectively; k is the 
viscosity rate constant; AE n and aE|< are viscous flow and cure activation 
energy, respectively; and n« and k* are two material constants. The 
viscosity is plotted as a function of time in Fig. 6. A linear dependence of the 
viscosity with respect to the logarithm of curing time has been commonly observed 
for many thermosetting resins. Values of parameters defined in Eq. (15) are 
tabulated below. 


Table I. Values of parameters defined in Eq. (18). 


n a = 3.78 x 10" 13 poises 
k = 5.5 x 10 9 min- 1 

CE 

AE = 28 Kcal/mole 
n 

AE k = 20 Kcal/mole 

R = 1.982 x 10“ 3 Kcal/mole °K 
T = 423°K 


By knowing resin chemoviscosity n(t) (Fig. 6) and fabric or fiber 
compression characteristics P G (h) (Fig. 3), it is possible to calculate the 
laminate thickness h(t) under a constant loading F by Eq. (10). The rate of 
decrease of the laminate thickness h(t) can be calculated by rearranging Eq. (9) 
so that 


dh _ 1 

it " nttT 


[ Kc z + f 0^ h3 l’ 


where and Q 2 are previously defined by Eq. (10). 


( 19 ) 
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Calculated values of h(t) are plotted in Fig. 7 for different values of 
C z Z 0 . As expected, the larger magnitudes of C 2 Z 0 result in thicker 
laminates. The rate of decrease in laminate thickness, dh/dt, is plotted for 
different values of C z Z 0 in Fig. 8. Similar to h(t), larger magnitudes of 
C z Z 0 are noted to result in lower values of dh/dt. Nearly linear relation- 
ships between ln(dh/dt) and t are observed for the range of values of 
C z Z 0 discussed herein. 

Values of the flow parameter Cp are also calculated and plotted against 
l-h/h 0 in Fig. 9 for different values of C z Z 0 . By comparing Figs. 5 and 
9 it can be seen that the present experiments have not reached the second plateau 
region (high Cp in Fig. 5) where fiber-fiber compression becomes dominant and 
the resin ceases to flow. It can be noted in Fig. 7 that the resin ceases to 
flow between 10 and 20 minutes after the load F = 1.39 x 10 8 dynes has been 
applied. This is due to the fast increase in chemoviscosity of the resin. 

In Fig. 10, l-h/h 0 is plotted against external load F for constant values 
of Cp according to Eq. (15). As F increases, l-h/h 0 approaches its 
limiting value of 1.0. A higher value of the flow parameter Cp implies lower 
viscosity resin, and l-h/h 0 approaches 1.0 at lower external loads, as 
expected. It should be emphasized that the quantity 1 - h(t)/h 0 , as defined 
in Fig. 2, is a measure of separation between two individual parallel prepreg 
layers at any instant during laminate processing, and is not a direct measurement 
of the change of thickness of the multi-layer laminate. The change of laminate 
thickness as a whole can, however, be calculated from (1 - h/h 0 ) by taking the 
fiber/fabric thicknesses into account. The curves shown in Fig. 10 are useful in 
establishing quality control of prepreg materials. For a given prepreg tape, 

Cp should be a constant if the load F is applied for a fixed period of time 
according to Eq. (11). The experiments to verify the curves in Fig. 10 are easy 
to implement in principle by means of a plastometer, and the validity of the 
assumption that C z Z 0 is a constant (instead of a function of h) throughout 
the tests can therefore be evaluated. 

A similar plot of l-h/h 0 versus F for different values of C z Z 0 is 
shown in Fig. 11. Because higher values of C z Z 0 result in lower values of 
dh/dt as discussed previously in regard to Fig. 8, the shifts of the curves to 
the higher F region on these plots as values of C z Z 0 become larger are 


expected. Figure 11 is constructed for a fixed value of Cp = 0.01, and is 
particularly useful as a complement to Fig. 10 in exploring validities of various 
assumptions used in the model. 

Resin losses in the vertical and horizontal directions can be calculated as 
well with this model. Fig. 12 shows the volumetric resin loss, Q z , in the 
vertical direction according to Eq. (4) for different values of C z Z 0 . It is 
obvious that the larger the pressure drop C z Z 0 is, the higher the resin loss 
in the z-di recti on becomes. The ratio of Q z to the total resin loss can be 
calculated by combining Eqs. (4) and (13). 

57 = - C Z K W dt/ / / dt ’ (20) 

where the second integral can be evaluated from Fig. 7. 
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IV. CONCLUSIONS 


A theory which describes the resin consolidation process during prepreg 
lamination has been formulated. It is based on the assumption that the system is 
composed of two types of flow: a horizontal squeezing flow between two nonporous 

parallel plates, and a vertical poiseuille type pressure flow through porous 
media. Flows in all directions originate from the squeezing action between two 
approaching fiber/fabric layers which are subjected to an external load F. A 
parametric study reveals that model predictions are in reasonable qualitative 
agreement with physical observations. Quantities such as resin loss during 
lamination, velocity and pressure profiles associated with resin flows, laminate 
consolidation speed and prepreg processibility, etc. , can all be calculated. The 
interrelationships between these quantities have been illustrated. Potential 
uses of the model in quality control of incoming prepreg materials have been 
discussed as well. Confirmations of the prepreg characterization methods, as 
suggested in the theory by experimental measurements, have been planned. The 
parallel -pi ate plastometer will be an appropriate device to verify the model. 
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Figure 1. Schematic diagram of multilayer prepreg. laminate 
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Figure 6. Chemoviscosity of a hypothetical B-stage 
epoxy resin at 150°C (423°K) 
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Figure 11. Changes of prepreg thickness vs. external load at different values 
of C Z 
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Figure 12. Volumetric resin flow in the vertical direction as a 
function of time at different values of C_Z_ 
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material are discussed. Experimental measurements needed to confirm the prepreg 
characterization methods suggested by the model have been planned for the future 
research. The parallel -pi ate plastometer will be an appropriate device to carry 
out the experiments. 
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